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Abstract 

A  one-dimensional,  two-phase,  transient  mathematical  model  was  developed  to  analyze  how  carbon  monoxide  poisoning  affects  the  performance 
of  a  PEM  fuel  cell.  This  work  examines  both  vapor  and  liquid  water  transport  inside  the  cell.  The  theoretical  results  indicate  that  a  higher  CO 
concentration  results  in  less  hydrogen  coverage  and  a  large  drop  in  the  time  to  reach  steady  state  tss.  The  slowing  of  the  reactions  at  both  the  anode 
and  the  cathode  reduce  the  saturation  of  liquid  water  in  the  catalytic  layers.  The  distribution  of  liquid  water  depends  more  strongly  on  the  CO 
concentration  than  on  dilution  of  hydrogen  in  the  MEA  of  the  fuel  cell.  Increasing  the  amount  of  pure  hydrogen  drastically  increases  tss  for  a  wide 
range  of  CO  contents.  At  a  relatively  low  CO  content,  a  long  tss  can  be  achieved  using  pure  hydrogen,  especially  at  high  cell  voltage,  promoting 
the  tolerance  for  CO  and  providing  the  desired  performance  of  the  fuel  cell. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  most  efficient  fuel  for  use  in  a  PEM  fuel  cell  is  pure 
hydrogen.  However,  this  is  difficult  to  store  and  has  a  high  cost  of 
production.  Reforming  from  hydrocarbons,  including  gasoline 
and  alcohol  is  the  most  extensively  used  technique  for  generating 
hydrogen  fuel  for  use  in  hydrogen  proton  exchange  membrane 
fuel  cells,  which  contain  45%  hydrogen,  10  ppm  CO,  15%  CO2 
and  1%  CH4  [1].  CO  at  a  concentration  even  as  low  as  5-10  ppm 
effectively  blocks  Pt  reaction  sites  for  the  electro-oxidation  of 
hydrogen  [2,3] .  The  affinity  between  CO  and  Pt  is  such  that  even 
at  this  low  concentration  of  CO,  the  consequent  loss  of  perfor¬ 
mance  is  severe  [4,5].  Lee  et  al.  [4]  investigated  the  tolerance  of 
PtSn/C,  Pt/C  and  PtRu/C  alloy  electrocatalysts  for  CO.  Depend¬ 
ing  on  the  nature  of  the  electrode  material,  the  CO  oxidation 
occurs  at  different  potentials.  The  improvement  of  CO  tolerance 
is  contributed  by  the  changes  in  the  thermodynamics  and  the 
CO  adsorption  process.  Murthy  et  al.  [5]  presented  the  steady 
state  and  transient  performance  of  a  fuel  cell  with  relatively  high 
concentrations  of  CO,  for  two  types  of  gas  diffusion  media. 
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Recently,  many  efforts  have  been  made  to  increase  the  tol¬ 
erance  of  the  PEM  fuel  cell  to  CO  [6-11].  Schmidt  et  al.  [6] 
and  Divisek  et  al.  [7]  presented  two  methods  for  improving  the 
cell  performance  using  H2/CO  as  a  fuel.  First,  the  use  of  Pt-Ru 
electrocatalysts  at  the  anode  can  considerably  enhance  the  tol¬ 
erance  to  CO.  Second,  the  addition  of  liquid  hydrogen  peroxide 
to  the  humidification  water  in  the  cell  leads  to  the  formation  of 
active  oxygen  by  the  decomposition  of  H2O2.  Complete  recov¬ 
ery  can  be  achieved  for  H2  100  ppm-1  CO.  Si  et  al.  studied  the 
CO  tolerance  of  the  Pt-Ru/C  catalyst  at  elevated  temperature  and 
atmospheric  pressure  in  a  PEM  fuel  cell  [8].  They  demonstrated 
that  the  anodic  polarization  fell  dramatically  as  the  temperature 
increased,  because  the  CO  coverage  was  lower.  Improving  the 
activity  of  the  catalyst  with  a  high  CO  tolerance  increases  the 
rates  of  reaction  of  CO  and  hydrogen.  Zhang  et  al.  [9]  found  that 
CO  poisoning  process  can  be  accelerated  at  high  anode  flow  rate. 
The  performance  decreases  substantially  as  the  anode  flow  rate 
increases,  because  CO  is  adsorbed  on  the  catalyst  site.  How¬ 
ever,  increasing  the  anode  overpotential  promotes  the  oxidation 
of  CO,  maintaining  desired  cell  performance.  Yu  et  al.  [10]  and 
Santiago  et  al.  [11]  adopted  the  Pt-Ru/C  electrocatalyst  for  dif¬ 
ferent  electrode  structures  to  improve  the  tolerance  of  PEMFC 
to  CO.  Water  sufficiently  activates  the  Ru  surface  to  promote  the 
oxidation  of  CO. 
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Nomenclature 

a  contact  area  of  Pt  catalyst  (cm2  cm-3) 

bm  ratio  of  forward  to  backward  of  hydrogen  adsorp¬ 
tion  (atm) 

bfco  ratio  of  forward  to  backward  of  CO  adsorption 
(atm) 

C  concentration  of  reactant  gas  (mol  cm-3) 

D  diffusion  coefficient  (cm2  s-1) 

i  current  density  (A  cm-2) 

k  conductivity  of  the  Nation  phase 

keco  CO  electro-oxidation  rate  constant  (A  cm-2) 
ke h  hydrogen  electro-oxidation  rate  constant 

(A  cm-2) 

kfco  CO  adsorption  rate  constant  (A  cm-2  atm- 1 ) 

km  hydrogen  adsorption  rate  constant 

(A  cm-2  atm-1) 
n  number  of  electrons 

P  total  pressure  (atm) 

R  universal  gas  constant  (J  mol-1  K-1) 

s  liquid  water  saturation 

T  temperature  (K) 

t  time  (s) 

X  molar  fraction 

v  distance  (pm) 

Greek  letters 

s  porosity 

0  ionic  potential  (V) 

y  stoichiometric  coefficient 

Oco  coverage  ratio  of  CO  on  Pt  catalyst  site 

On  coverage  ratio  of  hydrogen  on  Pt  catalyst  site 

§  molar  area  density  of  Pt  catalyst  sites  (C  cm-2) 

Superscripts 
CL  catalyst  layer 

in  inlet  at  catalyst  layer 

M  membrane 

Subscripts 

CO  carbon  monoxide 

H2  hydrogen 

O2  oxygen 

ss  steady  state 

wg  vapor  water 

wn  liquid  water  in  Nafion  phase 


In  theoretical  studies,  Springer  et  al.  [12]  derived  a  mathe¬ 
matical  model  to  describe  CO  poisoning  on  the  catalytic  sites. 
Chan  et  al.  [13]  combined  the  theoretical  models  developed 
by  Springer  et  al.  [12]  and  Bernardi  and  Verbrugge  [14,15] 
to  examine  the  CO  kinetics.  Bhatia  and  Wang  [16]  treated  the 
characteristics  of  the  anode  catalyst  layer  as  a  boundary  con¬ 
dition  in  analyzing  the  transient  CO  poisoning  behaviors  for 
various  levels  of  CO.  Baschuk  and  Li  [17]  developed  a  math¬ 


ematical  model  to  simulate  both  CO  poisoning  and  oxygen 
bleeding. 

Although  the  fuel  cell  performance  exhibited  transient  varia¬ 
tions  at  different  CO  concentrations,  the  actual  coverage  profile 
and  the  reactant  gas  distribution  over  the  catalyst  layer  are  not 
presented.  Even  the  transportation  of  liquid  water  and  its  dis¬ 
tribution  in  the  MEA  during  transient  CO  poisoning  are  not 
well  understood.  In  this  work,  a  transient  two-phase  mathemat¬ 
ical  model  of  the  poisoning  effect  of  PEM  fuel  cells  by  CO  is 
investigated.  Various  parameters  are  considered  to  promote  the 
tolerance  for  CO  and  analyze  the  effect  on  the  characteristic  time 
ts s,  and  thus  elucidate  transport  phenomena  inside  the  cell  and 
improve  its  performance. 

2.  Theoretical  model 

Fig.  1  depicts  the  presented  simulation  model  of  the  PEM 
fuel  cell,  including  the  anode  catalyst  layer,  the  membrane  and 
the  cathode  catalyst  layer.  The  transient  two-phase  theoretical 
CO  poisoning  behavior  is  investigated.  Table  1  presents  the 
governing  equations,  where  gaseous  hydrogen  concentration 
Ch?  ;  the  gaseous  oxygen  concentration  Co2 ;  the  gaseous  carbon 
monoxide  concentration  Cqo\  the  vapor  concentration  Cwg;  the 
concentration  of  liquid  water  in  the  Nafion  phase  Cwn;  the  satu¬ 
ration  of  liquid  water  s  and  the  ionic  potential  0.  The  governing 
equations  must  be  solved  simultaneous  for  the  dependent  vari¬ 
ables.  The  reactant  gas  and  vapor  water  transport  are  expressed 
as: 

Ni  =  -Di[sCL0  -s)]L5VCi  (1) 

where  £CL  is  the  total  porosity  and  s  is  the  liquid  water  saturation 
on  both  the  anode  and  the  cathode  sides.  The  water  transfer  rate 
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Table  1 


Governing  equations 


Variables 

ACL 

MEM 

CCL 

Ch2 

§  [(1  ~  V)£CLCh2]  =  — V  '  ^h2  ~  ^h2 

Cr2  =  0 

Cr2  =  0 

Co2 

C02  =0 

C02  =0 

f[(l-s)£CLCo2]  =  -V./Vo2- 

Ro2 

Cco 

|[(1  -  s)£clCco]  =  -V  •  Nco  -  Rco 

Cco  =  0 

Cco  =  0 

Cwg 

97  [(1  —  =  —V  •  Nw g  —  Rw 

CWg  =  0 

|[(l-^CLCwg]  =  -V-/Vwg- 

/?w 

GWn 

^  Ovn  =  0 

9Cwn  _  /-»  u2  n 

c  —  ^wn  v  ^wn 

^  Cw  n  =  0 

CL 

C  %  =  ;°  (  d£ )  O*2* + ■  p2) + Rw 

s  =  0 

CL 

£0  Pw  ds  _  Pw^w,0  /  d Pc\  (  V72 
Mw  dt  ~  Mw/xw  V  d,s’  )  V 

(An +  2 ^  Ro2  +  Rw 

+  (V  ■  s)2)  + 

0 

^effV20  —  (/?h2  +  ^co)  =  0 

&V2</>  =  0 

^effV20  -  Rq2  =  0 

between  the  gas  and  liquid  interfacial  phase  used  by  Lin  et  al. 
[19]  and  He  et  al.  [18]  is  used  herein: 


dOpp 

d  t 


=  ^fCO^CO^(l 


Oh  -  Oco )  -  ^fco^fco^co 


flw  =kc g°  (lRTS)yv(yvP-  rf)q 

CL 

+kv  (yv  p  -  pf)(  i  -  9)  (2) 

The  first  and  second  term  represents  the  condensation  and  evap¬ 
oration  rates,  respectively,  where  q  is  the  switching  function 
while  kc  and  kv  are  rate  constants.  The  transport  of  liquid  water 
in  the  membrane  is  driven  by  the  combined  effect  of  diffusion 
and  electro-osmotic  drag  [20] : 

9Cwn  i 

Nw  =  -Dwn—^  +  -nd  (3) 

ox  F 

in  which  Nw  is  the  flux  of  liquid  water,  Dwn  the  diffusion  coef¬ 
ficient  of  liquid  water,  i  the  flux  of  the  charges,  nd  the  electro- 
osmotic  drag  coefficient,  and  F  is  the  Faraday  constant.  The 
transport  of  liquid  water  in  the  catalytic  layer  is  described  by 
the  simplified  correlation  Kw(s)  =  Kw$s,  which  elucidates  the 
relationship  between  the  permeability  to  liquid  water  and  the 
saturation  of  liquid  water  [19]. 

In  the  anode  catalytic  layer,  the  distributions  of  the  current 
densities  of  hydrogen  and  carbon  monoxide  are: 


d /h2 

dx 


=  2<2&eH$H  h 


nu2  F{4>s  -  0 


2  RT 


d/co 

dx 


=  2akepp0pp  sin  h 


f  nppF((ps  -  0 
V  2  RT 


where  ke h  and  kQpp  are  the  hydrogen  and  CO  electro-oxidation 
rate  constants  [12],  <fis  the  electronic  potential,  0  the  electrolytic 
phase  potential  [21]  and  Uq  is  the  thermodynamic  equilibrium 
potential.  The  coverage  Oh  and  Opp  is  [12]: 


dOn 

dt 


=  ^fh^hFCI 


Oh  Oqo ) 


— 2fceH$H  sin  h 


f  nn2  F(</>s  ~  0 
V  2  RT 


—2kQpp0pp  exp 


(  ncoFjcps  -  0 

V  2  RT 


where  §  is  the  molar  area  density  of  the  catalyst  sites,  Xh  the 
molar  fraction  of  hydrogen,  Xco  the  CO  molar  fraction,  and  P  is 
the  total  pressure.  The  reaction  rates  of  hydrogen  and  CO  within 
the  anode  catalyst  layer  are: 


) H2  /  Kh2  \ 

dx  \  nn2  F  J 


Rco 


Yco 

nppF 


where  y  is  the  stoichiometric  coefficient  and  n  is  the  number 
of  electrons.  At  the  cathode  catalyst  layer,  the  reaction  rate  of 
oxygen  is  given  by: 


=  j=(l  ~  0co)(l  -  s)aio 
4  F 


x 


C02 

- exp 

_Q)2,ref 


^-f(0,  -  0-  cQ)y 


(10) 


where  Opo  is  the  mean  CO  coverage,  which  was  calculated  from 
the  anodic  CO  coverage  0pp.  The  factor  (1  —  s)  represents  the 
effect  of  liquid  water  saturation  in  the  catalyst  layers.  The  initial 
conditions  are  all  steady  state  without  zero  CO.  Table  2  lists 
all  of  the  boundary  conditions  used  in  this  simulation  model. 
Without  consideration  for  the  gas  diffusion  layer,  the  liquid  water 
saturation  is  assumed  to  be  zero  at  the  boundaries  (x  =  0  and  L). 


3.  Results  and  discussion 


Various  CO  concentrations  and  hydrogen  dilutions  are 
employed  to  simulate  a  wide  range  of  hydrogen  fuels  from  the 
reformer,  and  thus  elucidate  the  transient  behavior  associated 
with  poisoning.  The  feed  vapor  water  is  fully  saturated  at  80  °C 
at  both  the  anode  and  the  cathode  inlets.  The  reactant  gas  distri¬ 
bution,  the  coverage,  the  liquid  water  distribution  and  the  time  to 
reach  steady  state  are  investigated.  The  initial  condition  is  started 
in  the  absence  of  CO.  The  steady-state  condition  is  defined  as 
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Table  2 


Boundary  conditions 


Variables 

x  =  0 

ACL/MEM 

MEM/CCL 

x-L 

Ch2 

Ch2  = 

A'-in 

ch2 

/Vh2Iacl  =  0 

N/A 

N/A 

co2 

N/A 

N/A 

No2  1  ACL  =  0 

Co2  =  cs2 

Cco 

Cco  = 

A'dn 

^CO 

NcoIacl  =  0 

N/A 

N/A 

Cwg 

Cwg  = 

/^in 

^wg 

N wg  1  ACL  =  0 

N wg  1 CCL  =  0 

c  —  rin 

<-Wg  —  '-wg 

Cwn 

N/A 

k(Ca  Cwn)  =  A/w  |mem 

2f  “I-  k(Cc  Cwn)  =  /Vw  MEM 

N/A 

s  =  0 

VvIacl  =/Vw|mem 

/VwImEM  =/V\vlcCL 

s  =  0 

0  =  0 

kn,t  ff  V0  =  knV0 

kn,e ff  V0  =  kn  V0 

v0  =  o 

pertaining  when  the  relative  error  of  the  output  current  density 
reaches  10-3  (A  cm-2): 

{i(n)-i(n  +  1)}  <  10"3  (11) 

in  which  i(n)  is  the  current  density  at  the  nth  min.  Table  3  presents 
the  parameters  used  in  this  work. 

Fig.  2  plots  the  transient  evolution  of  the  hydrogen  cover¬ 
age  across  the  anode  catalyst  layer  with  10  ppm  CO,  100%  H2 
and  0.6  V.  The  simulation  begins  in  the  steady  state  without  CO 
poisoning.  Consequently,  the  hydrogen  coverage,  On,  decreases 
with  time  because  CO  is  adsorbed  on  the  Pt  catalyst,  reducing 
the  number  of  catalytic  sites  available  for  the  electro-oxidation 
of  hydrogen.  Fig.  3  shows  that  during  poisoning  CO  is  adsorbed 
on  the  catalytic  sites,  because  the  Pt  catalyst  has  a  strong  affin¬ 
ity  for  CO.  The  accumulation  of  CO  at  the  catalytic  sites  is 
sustained,  inhibiting  the  electro-oxidation  of  hydrogen.  Hence, 
the  hydrogen  and  CO  coverage  takes  100  min  to  reach  a  steady 
state. 

Fig.  4  shows  the  liquid  water  saturation  profiles  across  the 
anode  catalyst  layer  at  different  times,  with  100%  H2,  10  ppm 
CO  and  0.6  V.  The  effect  of  the  electro-osmotic  drag  is  pro¬ 
portional  to  the  current  density,  which  was  generated  by  the 
electro-chemical  reaction  of  hydrogen  and  oxygen.  The  cov¬ 
erage  of  hydrogen  falls  with  time,  so  the  current  density  was 
reduced,  weakening  the  effect  of  electro-osmotic  drag.  However, 

Table  3 

The  parameters  used  in  the  present  model 


Temperature,  T  (K) 

Total  pressure,  P  (atm) 

Diffusion  coefficient  of  hydrogen  in  gas 
phase,  Dho  (cm2  s-1) 

Diffusion  coefficient  of  oxygen  in  gas  phase, 
D0l  (cm2  s— 1 ) 

Diffusion  coefficient  of  vapor  in  gas  phase, 
Dw g  (cm2  s-1) 

Thickness  of  catalyst  layer,  <5cl  (pan) 
Thickness  of  membrane,  <5mem  (pan) 

Gas  porosity  in  catalyst  layer,  eCL 
Volumetric  fraction  of  Nation  in  membrane 
membrane,  £m 

Ionic  conductivity,  kn  (mho  cm-1) 
hydrogen  adsorption  rate  constant,  km 
(A  cm-2  atm-1) 

ratio  of  forward  to  backward  of  CO 
adsorption,  bfco  (atm) 


353 

1 

1.1028 

0.1775  x  (77273.15)1-823 

0.256  x  (77307. 15)2-334 

16 

50 

0.4 

0.4 

0.17 

100 

1.7  x  10-6 


Fig.  2.  The  transient  evolution  of  the  hydrogen  coverage  profile  across  the  anode 
catalyst  layer,  with  100%  H2,  10  ppm  CO,  0.6  V. 


Fig.  3.  The  transient  evolution  of  the  CO  coverage  profile  across  the  anode 
catalyst  layer,  with  100%  H2,  10  ppm  CO,  0.6  V. 
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Anode  Catalyst  Layer,  pm 


Fig.  4.  The  transient  evolution  of  the  liquid  water  saturation  profile  across  the 
anode  catalyst  layer,  with  100%  H2,  10  ppm  CO,  0.6  V. 

the  oxygen  reduction  reactions  are  also  suppressed,  reducing 
the  diffusion  of  water  from  the  cathode  to  the  anode.  Therefore, 
the  saturation  by  liquid  water  across  the  anode  catalyst  layer 
decreases  with  time.  During  poisoning,  the  diffusion  of  liquid 
water  from  cathode  to  anode  dominates  the  amount  of  liquid 
water  at  saturation  in  the  anode  catalyst  layer. 

Fig.  5  plots  the  unsteady  variations  of  the  liquid  water  dis¬ 
tribution  across  the  membrane  at  various  times.  The  gradient  of 
the  liquid  water  distribution  declines  with  time,  because,  as  the 
duration  of  poisoning  increases,  the  rate  of  the  reaction  decreases 
with  time  on  both  the  anode  and  the  cathode  sides.  The  effect  of 
electro-osmotic  drag  and  the  diffusion  of  liquid  water  from  the 
cathode  to  the  anode  are  also  weakened,  reducing  the  slope  of 


Fig.  5.  The  transient  evolution  of  the  water  content  profile  across  the  membrane, 
with  100%  H2,  10  ppm  CO,  0.6  V. 


Fig.  6.  The  transient  evolution  of  the  liquid  water  saturation  profile  across  the 
cathode  catalyst  layer,  with  100%  H2,  10  ppm  CO,  0.6  V. 

the  liquid  water  distribution  across  the  membrane.  Fig.  6  reveals 
that  the  amount  of  liquid  water  at  saturation  of  the  cathode  cata¬ 
lyst  layer  greatly  exceeds  that  on  the  anodic  side  of  the  catalyst. 
The  small  electro-osmotic  drag  and  the  generation  of  less  liquid 
water  cause  the  saturation  level  to  drop  with  time. 

Fig.  7  presents  the  ionic  potential  profile  across  the  ME  A 
at  various  times.  No  gas  fuel  is  consumed  in  the  membrane 
phase,  so  the  ionic  potential  distribution  plotted  in  Fig.  7  is  a 
straight  line.  In  the  catalytic  layers,  and  especially  in  the  cath¬ 
ode  catalyst  layer,  the  ionic  potential  distribution  is  nonlinear, 
because  reduction  reactions  consume  oxygen.  During  poison¬ 
ing,  the  ionic  potential  falls  with  time  reducing  the  output  current 
density. 


Fig.  7.  The  transient  evolution  of  the  ionic  potential  profile  across  the  MEA, 
with  100%  H2,  10  ppm  CO,  0.6  V. 
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Anode  Catalyst  Layer,  pm 

Fig.  8.  The  steady  state  hydrogen  coverage  profile  across  the  anode  catalyst 
layer  at  various  CO  concentrations  and  hydrogen  dilutions  at  0.6  V. 


Fig.  10.  The  steady  state  liquid  water  saturation  profile  across  the  anode  catalyst 
layer  at  various  CO  concentrations  and  hydrogen  dilutions  at  0.6  V. 


Fig.  8  depicts  the  steady  state  hydrogen  coverage  profile 
across  the  anode  catalyst  layer  at  various  CO  concentrations  and 
hydrogen  dilutions  at  0.6  V.  The  results  indicate  that  at  10  ppm 
CO,  tss  is  longer  than  at  100  ppm  CO,  for  both  40%  H2  and  100% 
H2 .  The  hydrogen  coverage  is  also  higher  at  lower  CO  concentra¬ 
tion,  because  increasing  the  CO  concentration  increases  the  rate 
of  adsorption  of  CO  onto  the  sites  of  the  catalyst.  Fig.  9  plots  the 
opposite  trend.  Hence,  hydrogen  dilution  significantly  affects 
the  hydrogen  coverage,  especially  at  low  ppm  CO.  During  poi¬ 
soning,  adding  pure  hydrogen  fuel  increase  more  reaction  sites 
for  hydrogen,  especially  at  low  CO  concentration.  As  shown  in 
Figs.  8  and  9,  when  an  anodic  inlet  flow  containing  an  H2/CO 
mixture  was  fed  into  the  fuel  cell,  CO  was  accumulated  on 
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Fig.  9.  The  steady  state  CO  coverage  profile  across  the  anode  catalyst  layer  at 
various  CO  concentrations  and  hydrogen  dilutions  at  0.6  V. 


the  catalytic  sites.  The  coverage  of  CO  is  much  higher  with 
100  ppm  CO,  results  higher  total  surface  coverage.  Nevertheless, 
reduces  effective  available  catalyst  sites  On  for  electro-oxidation 
of  hydrogen. 

Fig.  10  plots  the  effect  of  CO  concentration  and  hydrogen 
dilution  on  the  distribution  of  liquid  water  saturation  across  the 
anodic  catalytic  layer.  Fig.  10  reveals  that  when  cells  are  operated 
with  10  ppm  CO,  the  amount  of  liquid  water  saturation  greatly 
exceeds  that  when  100  ppm  CO  is  used,  at  both  40%  and  100% 
hydrogen.  The  presence  of  CO  in  the  anodic  inlet  flow  inhibits 
the  reduction  reactions  on  cathode  side  and  declines  the  diffusion 
of  water  from  the  cathode  to  anode.  The  distribution  of  liquid 
water  saturation  depends  more  strongly  on  the  CO  concentra¬ 
tion  than  on  dilution  of  hydrogen  in  the  anodic  catalytic  layer. 
Fig.  1 1  plots  the  liquid  water  distribution  across  the  membrane. 
At  10  ppm  CO,  100%  H2,  the  effect  of  the  electro-osmotic  drag  is 
strong  and  much  liquid  water  is  generated  at  the  cathode  catalyst 
layer,  causing  the  liquid  water  distribution  to  have  a  large  slop. 
Consequently,  the  CO  concentration  significantly  influences  the 
distribution  of  liquid  water  across  the  membrane.  The  dilution 
of  hydrogen  by  the  fuel  reduces  the  gradient  of  the  liquid  water 
distribution  because  the  electro-osmotic  drag  is  small  and  less 
water  is  generated  in  the  cathode  catalyst  layer.  An  increase  in 
the  CO  concentration  effectively  reduces  the  gradient  of  liquid 
water  distribution  across  the  membrane. 

Fig.  12  presents  the  profile  of  liquid  water  saturation  across 
the  cathode  catalyst  layer.  This  figure  shows  that  an  increase 
in  CO  concentration  markedly  affects  the  liquid  water  satu¬ 
ration  profile  because  the  oxygen  reduction  reaction  drops  in 
the  cathode  catalyst  layer.  As  the  CO  level  increases,  or  the 
hydrogen  becomes  more  dilute,  the  saturation  profile  becomes 
lower  shown  in  Fig.  12,  the  CO  content  strongly  affects  the  liq¬ 
uid  saturation  at  the  cathode  catalyst  layer.  Fig.  13  plots  the 
ionic  potential  at  various  CO  levels  and  hydrogen  dilutions.  As 
discussed  above,  the  loss  of  ionic  potential  reduces  the  output 
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Fig.  11.  The  steady  state  water  content  profile  across  the  membrane  at  various 
CO  concentrations  and  hydrogen  dilutions  at  0.6  V. 

current  density.  Therefore,  increasing  the  CO  level  or  diluting 
the  hydrogen  seriously  reduces  the  loss  of  ionic  potential. 

Fig.  14  compares  the  presented  simulation  results  with  exper¬ 
imental  data  reported  by  Bhatia  and  Wang  [16].  There  are  four 
gas  compositions  fed  into  the  anode.  The  simulation  results  show 
that  the  cell  performance  and  durability  depends  strongly  on 
the  dilution  of  hydrogen  and  CO  concentration,  considerably 
reducing  the  current  density  and  the  time  tss  required  to  reach  a 
steady  state.  Increasing  hydrogen  dilution  and  CO  concentration 
degrades  the  performance  of  the  cell  and  reduces  the  time  fss.  At 
a  low  level  of  CO  (10  ppm  CO),  the  cell  performance  depends 
more  strongly  on  the  dilution  of  hydrogen,  considerably  reduc¬ 
ing  the  current  density.  The  predicted  CO  poisoning  results  agree 


Fig.  12.  The  steady  state  liquid  water  saturation  profile  across  the  cathode  cat¬ 
alyst  layer  at  various  CO  concentrations  and  hydrogen  dilutions  at  0.6  V. 


Fig.  13.  The  steady  state  ionic  potential  profile  across  the  MEA  at  various  CO 
concentrations  and  hydrogen  dilutions  at  0.6  V. 

closely  with  the  experimental  values,  except  from  the  case  for 
10  ppm,  40%  H2.  The  discrepancy  may  be  attributed  without 
consideration  for  the  gas  diffusion  layer.  The  assumption  may 
probably  influence  the  distribution  of  liquid  water  in  the  catalyst 
layers.  This  may  explain  the  discrepancy  between  the  predicted 
simulation  results  and  experiments. 

Fig.  15  plots  the  influence  of  CO  concentration  on  the  time 
to  reach  the  steady  state  for  two  hydrogen  dilutions  and  cell 
voltages.  The  theoretical  results  indicate  that  a  higher  CO  con¬ 
centration  results  in  large  drop  in  the  time  to  reach  steady  state 
tss.  The  dropping  rate  of  tss  becomes  small  at  a  large  CO  con¬ 
tent.  Hydrogen  dilution  substantially  increases  time  tss  under  a 
wide  range  of  CO  content.  At  high  cell  voltages,  a  large  time 
tss  is  obtained,  especially  at  small  CO  content,  because  only  a 
small  amount  of  hydrogen  is  used  at  low  current  density,  so  its 


Fig.  14.  The  present  simulation  results  compared  with  experimental  data  at 
0.6  V. 
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Fig.  15.  The  influence  of  CO  concentration  on  the  time  to  reach  steady  state  for 
two  hydrogen  dilutions  and  cell  voltages. 

effect  on  the  CO  adsorption  rate  in  the  anode  catalyst  layer  is 
weak.  However,  the  cell  voltage  only  weakly  influences  the  time 
ts s  at  high  CO  content.  In  this  study,  increasing  the  amount  of 
pure  hydrogen  drastically  increases  tss  for  a  wide  range  of  CO 
contents. 

The  results  for  pure  and  40%  hydrogen  at  two  CO  contents 
are  plotted  for  comparison  in  Fig.  16.  A  higher  CO  content  cor¬ 
responds  to  lower  tss.  At  100  ppm  CO,  the  cell  voltage  does  not 
clearly  affect  tss.  A  large  time  tss  can  be  achieved  using  pure 
hydrogen.  At  10  ppm  CO,  a  small  difference  on  tss  is  observed 
at  low  cell  voltage  for  pure  and  40%  hydrogen.  Thereafter,  the 
time  fss  increases  markedly  as  the  cell  voltage  increases  for  pure 
and  40%  hydrogen  at  10  ppm  CO.  Fig.  16  reveals  that  a  long 


Fig.  16.  The  effect  of  cell  voltage  on  the  time  to  reach  steady  state  for  two 
hydrogen  dilutions  and  CO  concentrations. 


4S  can  be  achieved  containing  pure  hydrogen  at  10  ppm  CO, 
especially  at  high  cell  voltage. 

4.  Conclusions 

This  study  investigated  the  transient  nature  of  poisoning  from 
carbon  monoxide  across  the  MEA  of  the  PEMFC  is  investigated. 
Hydrogen  dilution  and  CO  content  drastically  affect  the  perfor¬ 
mance  and  durability  of  the  fuel  cell.  When  an  anodic  inlet  flow 
containing  an  H2/CO  mixture  was  fed  into  the  fuel  cell,  CO  is 
accumulated  on  the  catalytic  sites.  This  is  because  the  Pt  catalyst 
has  a  strong  affinity  for  CO,  inhibiting  the  electro- oxidation  of 
hydrogen.  The  hydrogen  coverage  and  liquid  water  saturation 
declines  as  the  reaction  proceeded.  The  gradient  of  the  liquid 
water  distribution  across  the  membrane  and  the  ionic  potential 
also  fell  with  time.  The  distribution  of  liquid  water  depends  more 
strongly  on  the  CO  concentration  than  on  dilution  of  hydrogen 
in  the  MEA  of  the  fuel  cell.  The  theoretical  results  indicate  that 
a  higher  CO  concentration  results  in  large  drop  in  the  time  to 
reach  steady  state  4S.  In  this  study,  increasing  the  amount  of  pure 
hydrogen  drastically  increases  tss  for  a  wide  range  of  CO  con¬ 
tents.  At  100  ppm  CO,  the  cell  voltage  does  not  clearly  affect 
4S .  A  large  time  tss  can  be  achieved  using  pure  hydrogen.  At 
low  CO  content,  the  influence  of  hydrogen  dilution  on  the  time 
4s  is  weak  at  low  cell  voltage.  Thereafter,  the  time  tss  increases 
markedly  as  the  cell  voltage  increases  at  low  CO  content.  When 
the  anode  inlet  flow  contains  pure  hydrogen,  a  long  tss  results, 
especially  at  low  ppm  of  CO  and  high  cell  voltage,  providing 
the  desired  performance  and  durability  of  the  fuel  cell. 
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